function dx=fun(t,x)
f1=-8/3*x(1)+x(2)*x(3);
f2=-10*x(2)+10*x(3);
f3=-x(1)*x(2)+28*x(2)-x(3);
dx=[f1;f2;f3];